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1 Attorney Docket No. 84280 

2 

3 A METHOD TO ESTIMATE THE MECHANICAL PROPERTIES OF A SOLID 

4 MATERIAL SUBJECTED TO ISONIFICATION 

5 

6 STATEMENT OF GOVERNMENT INTEREST 

7 The invention described herein may be manufactured and used 

8 by or for the Government of the United States of America for 

9 governmental purposes without the payment of any royalties 
10 thereon or therefor. 

11 

12 BACKGROUND OF THE INVENTION 

13 (1) Field of the Invention 

14 The present invention relates to a method for measuring 

15 mechanical characteristics of materials. More particularly, 

16 this invention provides a method which uses transfer functions 

17 obtained by insonifying the material at different angles. Once 

18 obtained, the transfer functions are manipulated to yield closed 

19 form values of dilatational and shear wavespeeds. The wavespeeds 

20 are combined to determine complex Lame constants, complex 

21 Young's modulus, complex shear modulus, and complex Poisson' s 

22 ratio for the material. 
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1 (2) Description of the Prior Art 

2 Measuring the mechanical properties of slab-shaped (i.e., 

3 plates) materials are important in that these parameters 

4 significantly contribute to the static and dynamic response of 

5 structures built with such materials. Resonant techniques have 

6 been used to identify and measure longitudinal properties for 

7 many years (See D.M. Norris, Jr., and W.C. Young, "Complex 

8 Modulus Measurements by Longitudinal Vibration Testing, " 

9 Experimental Mechanics, Volume 10, 1970, pp. 93-96; W.M. 

10 Madigosky and G.F. Lee, "Improved Resonance Technique for. 

11 Materials Characterization," Journal of the Acoustical Society 

12 of America, Volume 73, Number 4, 1983, pp. 1374-1377; S.L. 

13 Garrett, "Resonant Acoustic Determination of Elastic Moduli," 

14 Journal of the Acoustical Society of America, Volume 88, Number 

15 1, 1990, pp. 210-220; G.F. Lee and B. Hartmann, U.S. Patent 

16 Number 5,363,701; G.W. Rhodes, A, Migliori, and R.D. Dixon, 

17 U.S. Patent Number 5,495,763; and R.F. Gibson and E.O. Ayorinde, 

18 U.S. Patent Number 5,533,399). 

19 These methods are based on comparing the measured 

20 eigenvalues of a structure to predicted eigenvalues from a model 

21 of the same structure. The model of the structure must have 

22 well-defined (typically closed form) eigenvalues for these 

23 methods to work. Additionally, resonant techniques only allow 

24 measurements at resonant frequencies. Most of these methods 
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1 typically do not measure shear wavespeeds (or modulus) and do 

2 not have the ability to estimate Poisson< s ratio. 

3 Comparison of analytical models to measured frequency 

4 response functions is another method used to estimate stiffness 

5 and loss parameters of a structure (See B.J. Dobson, "A 

6 Straight-Line Technique for Extracting Modal Properties From 

7 Frequency Response Data," Mechanical Systems and Signal 

8 Processing, Volume 1, 1987, pp. 29-40; T. Pritz, "Transfer 

9 Function Method for Investigating the Complex Modulus of 

10 Acoustic Materials: Rod-Like Specimen," Journal of Sound and 

11 Vibration, Volume 81, 1982, pp. 359-376; W.M. Madigosky and G.F. 

12 Lee, U.S. Patent Number 4,352,292; and W.M. Madigosky and G.F. 

13 Lee, U.S. Patent Number 4,418,573). When the analytical model 

14 agrees with one or more frequency response functions, the 

15 parameters used to calculate the analytical model are considered 

16 accurate. If the analytical model is formulated using a 

17 numerical method, a comparison of the model to the data can be 

18 difficult due to dispersion properties of the materials. 

19 Another method to measure stiffness and loss is to deform 

20 the material and measure the resistance to the indentation (See 

21 W.M. Madigosky, U.S. Patent Number 5,365,457). However, this 

22 method can physically damage the specimen if the deformation 

23 causes the sample to enter the plastic region of deformation. 
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1 Others methods have used inspnif ication as a means to 

2 determine defects in composite laminate materials (See D.E. 

3 Chimenti and Y. Bar-Cohen, U.S. Patent Number 4,674,334). 

4 However, these methods do not measure material properties. 

5 A method does exist to measure shear wave velocity and 

6 Poisson's ratio in the earth using boreholes and seismic 

7 receivers (See J.D. Ingram and O.Y. Liu, U.S. Patent Number 

8 4,633,449). However, this method needs a large volume of 

9 material and is not applicable to slab-shaped samples. 

10 Additionally, it needs a borehole in the volume at some 

11 location. 

12 in view of the above, there is a need for a method to 

13 measure complex frequency-dependent dilatational and shear 

14 wavespeeds of materials subject to insonif ication. Once the 

15 wavespeeds are identified, the complex frequency-dependent 

16 Young's and shear moduli and complex frequency-dependent 

17 Poisson's ratio can also be measured (or estimated). 
18 

19 SUMMARY OF THE INVENTION 

20 Accordingly, it is a general purpose and primary object of 

21 the present invention to provide a method to measure (or 

22 estimate) the complex frequency-dependent dilatational and shear 

23 wavespeeds of a slab of material subjected to insonif ication. 
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1 It is a further object of the present invention to provide 

2 a method to measure (or estimate) the shear modulus of a slab of 

3 material subjected to insonif ication . 

4 It is a still further object of the present invention to 

5 provide a method to measure (or estimate) the Young's modulus of 

6 a slab of material subjected to insonif ication. 

7 It is a still further object of the present invention to 

8 provide a method to measure (or estimate) the complex frequency- 

9 dependent Poisson' s ratio of a slab of material subjected to 
10 insonif ication. 

n to attain the objects described, there is provided a method 

12 which uses three transfer functions that are obtained by 

13 insonif ying the material at different angles. Once this is 

14 accomplished, the transfer functions are manipulated with an 

15 inverse method to yield closed form values of dilatational and 

16 shear wavespeeds at arty given test frequency. The wavespeeds are 

17 combined to determine complex Lame constants, complex Young's 

18 modulus, complex shear modulus, and complex Poisson' s ratio. 
19 

20 BRIEF DESCRIPTION OF THE DRESSINGS 

21 A more complete understanding of the invention and many of 

22 the attendant advantages thereto will be readily appreciated as 

23 the same becomes better understood by reference to the following 
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1 detailed description when considered in conjunction with the 

2 accompanying drawings wherein: 

3 FIG. 1 depicts a test setup to insonify and gather 

4 measurements for a specimen of material; 

5 FIG. 2 depicts the coordinate system of the test setup of 

6 FIG. 1; 

7 FIG. 3 is a plotted graph depicting the measurable transfer 

8 function of magnitude versus the frequency; 

9 FIG. 4 is a plotted graph depicting the measurable transfer 

10 function of phase angle versus the frequency; 

11 FIG. 5 is a plotted graph depicting the measurable transfer 

12 function "s" versus the frequency; 

13 FIG. 6 is a plotted graph depicting the real component of 

14 the actual and estimated dilatational wavespeed versus the 

15 frequency; 

16 FIG. 7 is a plotted graph depicting the imaginary component 

17 of the actual and estimated dilatational wavespeed versus the 

18 frequency; 

19 FIG. 8 is a plotted graph depicting the surface defined in 

20 equation (64) of the description versus the real and imaginary 

21 components of fi 2 at 1800 Hz with the magnitude depicted as a 

22 gray scale image; 

23 FIG. 9 is a plotted graph depicting the contour of the 

24 surface versus both the real and imaginary parts of p 2 ; 
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1 FIG. 10 is a plotted graph depicting the actual shear 

2 wavespeed and the estimated shear wavespeed versus the frequency 

3 with the real component; 

4 FIG. 11 is a plotted graph depicting the actual shear 

5 wavespeed and the estimated shear wavespeed versus the frequency 

6 with the imaginary component; 

7 FIG. 12 is a plotted graph depicting the actual shear 

8 modulus and the estimated shear modulus versus the frequency 

9 with the real component; 

10 FIG. 13 is a plotted graph depicting the actual shear 

11 modulus and the estimated shear modulus versus the frequency 

12 with the imaginary component; 

13 FIG. 14 is a plotted graph depicting the actual Young's 

14 - modulus and the estimated Young's modulus versus the frequency 

15 with the real component; 

16 FIG. 15 is a plotted graph depicting the actual Young's 

17 modulus and the estimated Young's modulus versus the frequency 

18 with the imaginary component; and 

19 FIG. 16 is a plotted graph depicting the actual Poisson's 

20 ratio and the estimated Poisson's ratio versus the frequency. 
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1 DESCRIPTION OF THE PREFERRED EMBODIMENT 

2 Referring now to the drawings wherein like numerals refer 

3 to like elements throughout the several views, one sees that 

4 FIG. 1 depicts the isonif ication of a slab-shaped test specimen 

5 10 by a speaker (or projector) 12. Insonif ication consists of 

6 loading the specimen 10 on one entire side with an acoustic wave 

7 originating at the speaker 12. The speaker 12 is located at a 

8 sufficient distance from the specimen 10 that the acoustic wave 

9 is nearly a plane wave by the time it contacts the specimen. 

10 The insonif ication is usually done at multiple frequencies and 

11 multiple angles. 

12 For the method presented, a frequency sweep (swept sine) is 

13 conducted at three different insonif ication angles. The 

14 transfer function data is collected with either accelerometers 

15 16, 18 on both sides which record accelerations, or laser 

16 velocimeters 20, 22 shining on both sides which record 

17 velocities. In the swept sine mode, the transfer functions of 

18 acceleration divided by acceleration or velocity divided by 

19 velocity are both equal to displacement divided by displacement. 

20 The time domain data are Fourier transformed into the frequency 

21 domain and then recorded as complex transfer functions, 

22 typically using a spectrum analyzer (not shown) . 
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I The motion of the specimen 10 is governed by the equation 

2 

3 ^V 2 nD + (A + //)VVo l m = p— 7- , U) 

a 

4 

2 

5 where A and are the complex Lame constants (N/m ),/>is the 

6 density (kg/m 3 ),o denotes a vector dot product; u is the 

7 Cartesian coordinate displacement vector of the material and d 

8 is the partial differential. 

9 The coordinate system of the test configuration is shown in 
10 FIG 2. Note that using this orientation results in b = 0 and a 

II having a value less than zero. The thickness h of the specimen 

12 10 is a positive value. Equation (1) is manipulated by writing 

13 the displacement vector u as 



14 



15 m = 



u x (x,y,z,t) 
u y {x,y,z,t) 
u z (x,y,z,t) 



(2) 



16 where x is the location along the specimen 10, y is. the location 

17 into the specimen 10, and z is the location normal to the 

18 specimen 10 and t is time (s) . The symbol V is the gradient 

19 vector differential operator written in three-dimensional 

20 Cartesian coordinates as 



21 
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2 

3 with i x denoting the unit vector in the x-direction, i y denoting 

4 the unit vector in the y-direction, and i z denoting the unit 

5 vector in the z-direction; V 2 is the three-dimensional Laplace 

6 operator operating on vector u as 



V 2 u = V 2 « x /, + V 2 m i + V 2 « 2 / 2 ' 



10 with V 2 operating on scalar u as 



(4) 



11 



12 



<? 2 « w d 2 u w d 1 *^ 



13 and the term V»u is called the divergence and is equal to 

14 

„ d*. d* y d* z ICS 
. V»u = — - + — - + . (6) 

16 The displacement vector u is written as 

17 

18 u = V^+Vx(? , (7) 
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1 where ^ is a dilatational scalar potential, x denotes a vector 

2 cross product, and ^ is an equivoluminal vector potential 

3 expressed as 



y/ x (x,y,z,t) 
y/ y (x,y,z,t) 



(8) 



The problem is formulated as a two-dimensional system, thus v=0, 
u y (x,y,z,t) = O r and d(-)/dy=0. Expanding equation (7) and breaking 
the displacement vector into its individual nonzero terms yields 



6 
7 
8 
9 
10 

11 

12 

13 and 

14 



u x (x,z,t) = 



d<fi(x,z,t) dy y {x,z,i) 



ex 



dz 



(9) 



Equations (9) and (10) are next inserted into equation (1), 
which results in 



15 

16 
17 

18 
19 

20 and 



(10) 



c 2 d V 2 0(x,z,t) = 



3 2 <f>(x,z,t) 



(11) 
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2 c f 2 VV,(x,z,0 = ^2 (12) 

3 

4 where equation (11) corresponds to the dilatational component 

5 and equation (12) corresponds to the shear component of the 

6 displacement field. Correspondingly, the constants c d and c s are 

7 the complex dilatational and shear wave speeds, respectively, 

8 and are determined by 



= i p 



c d = 



11 

12 and 
13 



15 

16 The relationship of the Lame constants to the Young's and shear 

17 moduli is shown as 
18 



19 



A = - (15) 

(l + v)(\-2u) 
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l and 

2 

„ = G = _J± — (16) 

3 M 2(l + y) ' 

4 where £ is the complex Young's modulus (N/m ), 6 is the complex 

5 shearmodulus (N/m 2 ), and u is the Poisson's ratio of the material 

6 (dimensionless) . 

7 The conditions of infinite length and steady-state response 

8 are now imposed, allowing the scalar and vector potential to be 

9 written as 



#(x, z,t) = <D(z) exp(iAx) exp(ifirf) (17) 



10 

11 

12 

13 ' and 

15 y/ y (x, z, t) = »F(z) exp(ifo) exp(ioO (18) 

16 

17 where i is the square root of -1, 6) is frequency (rad/s) , and k 

18 is wavenumber with respect to the x axis (rad/m),4> is the 

19 amplitude of the scalar potential 0 as a function of z and V is 

20 the amplitude of the vector potential as a function of z. 

21 Inserting equation (17) into equation (11) yields 
22 
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1 



^- + a <D(z) = 0 ' 



(19) 



2 where a is the modified dilatational wavenumber and d is the 

3 differential operator and where 



5 a = 4k]-k 2 , (20) 

6 

7 with k d is the actual dilatational wave number and 
8 

9 K=—. (2D 

10 

11 Inserting equation (18) into equation (12) yields 
12 

13 f^£> +/ ^( z ) = 0 > «22) 

dz 

14 

15 where 
16 

17 P = 4k]-k' 1 , (23) 

18 
19 

20 
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l with jS as the modified shear wavenumber and 



k =— . {24) 



4 

5 The solution to equation (19) is 

6 

7 O(z) = A{k, <y)exp(iaz) + B(k, <y)exp(-iaz) , (25) 

8 and the solution to equation (22) is 

9 4>(z) = C(Jfc, o)exp(iyfe) + D(k, <y)exp(-i/fe) , (26) 

10 

11 where i4(Jt,<»), B(k,0)), C(k,a>) , and are wave response 

12 coefficients that are determined below. The displacements can 

13 now be written as functions of the unknown constants using the 

14 expressions in equations (9) and (10) . They are 
15 

u z (x, z, 0 = U z (k, z, 6>)exp(iAx)exp(iatf ) 

16 r 1 > ( 27 > 

= { ia[A(k,co)exp(iaz) - B(k,G>)exp(-iaz)\+ 
\k[C{k, a) exp(iy?z) + D(k, <y)exp(-iy&)] )exp(ifcc)exp(ifijr) 

18 with [7 Z as the amplitude of displacement in the "z" direction and 
19 
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u x (x, z, t) = U x (k, z, <y)exp(ifa)exp(i<yO 

1 i < 28) 
= { ik[A(k,G))exp(iaz) + B(k,a>)^p(-iaz)\- 

iy5[C(it,«)exp(i>9z) - D(k, o)exp(-i^)] }exp(ifo)exp(ifflr) 

2 with U x as the amplitude of displacement in the "x" direction. 

t 

3 The normal stress at the top of the plate (z = b) is equal to 

4 opposite the pressure load created by the projector. This 

5 expression is 
6 

du,(x,bJ) .tfaAx,b,t) , . ,, Q v 

8 

9 and the tangential stress at the top of the plate b is zero and 
10 this equation is written as 
11 



12 T a (x,b,t) = fi 



~ai x (x,b,t) [ ai z (x,b,t) 
dz ' dx. 



= 0. (30) 



13 

14 The normal stress the bottom of the plate (z = a) is equal to 

15 zero. This expression is 
16 

„ ~ ^duAx,a,t) .duAx,a,i) 
17 Ta (x,a,t) = (A + 2 M ) ^ = 0 > (31) 

18 
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and the tangential stress at the bottom of the plate is zero and 
this equation is written as 



T a (x,a,t) = M 




(32) 



The applied load in equation (29) is an acoustic pressure 
and is modeled as a function at definite wavenumber and 
frequency as 

p 0 (x,z,t) = P 0 (»)exp(iAx)exp(ietf) , (33) 

with P being the amplitude and where the wavenumber k is found 
using 



where Cf is the compressional wavespeed of air (m/s) andflis the 

angle of incidence of the projector with the z axis (rad) . 

Assembling equations (1) - (34) and letting b = 0 yields 
the "A" matrix, x vector, and b vector in a four-by-four system 
of linear equations that model the system written in matrix 
form. They are 



k = — sin(0) , 

c f 



(34) 



17 



Ax = b ; (A 4x4, x 4xl f b 4x1) 



where the entries of equation (35) are 



A2 ~ Al » 



A l3 = 2kfifi , 



A4 ~ A3 > 



= -2/jka , 



A2 — Ai > 



^24 - ^23 > 



= 4,exp(iaa) , 
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2 
3 



A 32 =A n exp(-iaa) , (45) 



A„ = -A n exp(i/fc) , (4 6) 



4 4 =4 3 exp(-iy&*) , (47) 



^ 41 =^,exp(iaa), (48) 



A n =-A 2X exp(-\aa) , (4 9) 



5 
6 
7 
8 
9 

10 N A A3 =A 23 exp(ifia) , (50) 

11 

12 ^ = A v exp(-i>S3) , (51) 

13 

14 fc u =-P 0 (a>), (52) 

15 
16 
17 

18 *3»=0, (54) 

19 and 

20 

21 b M =0. (55) 



ft 2t =0, (53) 
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2 Using equations (35) - (55) the solution to the constants A(k,co) , 

3 B(k,6)), C(k 9 a)), and D{K(o) can be calculated at each specific 

4 wavenumber and frequency using 

5 s = A" l Ib . (56) 

6 

7 Once these are known, the transfer function T between the 

8 wall motion in the z direction at z = a and the wall motion in 

9 the z direction at z = b is now written in closed form notation 
10 using equations (27) and (56). The resulting expression is 

11 



12 



_ UA^aM _ *c*Pk 2 sin(<y/Q + Qg z -k z Y sin(^) (5?) 

( ,6)) U 2 (k 9 b,a>) 4^ 2 sm(a/0cosW 



13 

14 The first step is to solve for the response at zero 

15 wavenumber, or what is typically referred to as broadside 

16 response, to determine the dilatational wavespeed. At zero 

17 wavenumber, the angle between the direction of propagation of 

18 the insonif ication energy and the z axis is zero. The response 

19 of the structure to broadside energy is comprised entirely of 

20 dilatational waves, i.e., no shear waves are excited at zero 

21 wavenumber. Furthermore, the transfer function given in 

22 equation (57) reduces to 
23 
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T ^=—r^= T ^T> (58) 



3 where T x (or R^) is the measurement data from the experiment with 

4 a insonification angle of zero and is typically a frequency- 

5 dependent complex number and the subscript 1 denotes the first 

6 experiment. Equation (58) can be expanded into real and 

7 imaginary parts and solved, resulting in a value for a x at every 

8 frequency in' which a measurement is made. The solution to the 
9- real part Re of a x is 
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11 Re(a,) = 



14 

15 
16 



— Arc cosO ) + — - n even 
2h 2h 

(59) 

— Arc cosf-s) + n odd, 

2h 2h 



12 

13 where 



= [ReW)] 2 +[InW] 2 -VfRe^)] 2 +[Im(/? 1 )] 2 } 2 -^[Re(^)] 2 >2[Im(/? 1 )] 2 -l} , (60) 



17 and n is a non-negative integer and the capital A denotes the 

18 principal value of the inverse cosine function. The value of n 

19 is determined from the function s, which is a periodically 

20 varying cosine function with respect to frequency. At zero 
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1 frequency, nis 0. Every time s cycles through ;r radians (180 

2 degrees), n is increased by 1. When the solution to the real 

3 part of «! is found, the solution to the imaginary part Im of a, 

4 is then written as 

5 

6 h l °Z<\ cos[Re(ai)h] s in[Re(a,)A]J 



7 



8 The real and imaginary parts of a, from equations (59) and (61) 

9 respectively are combined to yield the complex wavenumber. 

10 Because this measurement is made at zero wavenumber (k = 0) , this 

11 is equal to the dilatational wavenumber. Thus, the dilatational 

12 wavespeed is equal to 
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14 C d = 



<» (62) 

[Re(a,) + ilm(a,)] 



15 

16 To solve for the shear wavespeed, the specimen must be excited 

17 at a nonzero wavenumber. This is done next. 

18 The next step is to solve for the response at nonzero 

19 wavenumber to determine the shear wavespeed. At nonzero 

20 wavenumber, the transfer function is given in equation (57). 

21 For this nonzero angle of insonif ication, this can be expressed 

22 as 



22 



0 



m 4a 2 fl*, 2 sm(a,h)HPl -*j) 2 MM = T 2 =— , (63) 

I l \K,G>) - 4a ^ k > sin(a2 / l)cos( ^/ l) + (> 9 2 2 -kl? cos(a 2 h)sHM R> 

2 

3 where T 2 (or R 2 ) is the measurement data from the experiment at 

4 nonzero insonif ication angle and is typically a frequency- 

5 dependent complex number and the subscript 2 denotes the second 

6 angle or experiment. It is noted that a 2 in equation (63) is 

7 different from a x in equation (58). This difference is based on 

8 a k 2 term shown in equation (20) where the wavenumber a is 

9 defined. Due to the complexity of equation (63), there is no 

10 simple method to rewrite the equation as a function f of P 2 > tne 

II variable that is to be estimated. Equation (63) can be 
12 rewritten as 

13 

14 /(/? 2 ) = 0 = 4a 2 /? 2 fc 2 2 sm(a 2 /^ > ( 64 > 

15 

16 where the problem now becomes finding the zeros of the right 

17 hand side of equation (64), or in the presence of actual data 

18 that contains noise, finding the relative minima of the right 

19 hand side of equation (64) and determining which relative 

20 minimum corresponds to shear wave propagation and which relative 

21 minima are extraneous. Because equation (64) has a number of 

22 relative minima, zero finding algorithms are not applied to this 

23 



1 function, as they typically do not find all of the minima 

2 locations. The best method to find all of the minima locations 

3 is by plotting the absolute value of the right hand side of 

4 equation (64) as a surface with the real part of p 2 on one axis 

5 and the imaginary part of p 2 .on the other axis. The value a 2 is 

6 determined using 



8 a 2 =Jkf^ = yl<z?-k 2 2 , (65) 

9 

10 so that equation (64) is a function of only /? 2 • Once this 

11 function is plotted, the minima can be easily identified and the 

12 corresponding value of (Jh) m at the location of the minima can be 

13 determined by examination of the minimum location, sometimes 

14 referred to as the grid method. The shear wave speed (s) are 

15 then determined using 
16 

17 



ou,*V(A)i+*2 (66) 



18 

19 and 
20 



21 (c,) m =— < 67 > 



22 
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1 where the subscript m denotes each minima value that was found 

2 by inspecting the surface formed from equation (64) . The 

3 determination of the correct index of m that corresponds to 

4 shear wave propagation is done below. 

5 The material properties such as Young's modulus and other 

6 material properties can be determined from the wavespeeds. The 

7 Lame constants are calculated with equations (13) and (14) 

8 written as 
9 

10 M m =P(cX (68) 

11 

12 and 
13 

14 X m =pc 2 d -2p{c s ) 2 m . (69) 

15 

16 To determine the correct index m that corresponds to the actual 

17 wave propagation rather than an extraneous solution, a third set 

18 of measurements are made at a nonzero incidence angle that is 

19 not equal to the angle used previously. The model in equation 

20 (63) is calculated from the estimated material properties and a 

21 residual value is determined using the third set of 

22 measurements. Each in indexed residual at a specific frequency 

23 is defined as 
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4a 3 (>g3)^3 2 sin(«3/») + [(>g3^ zMl^iMS J_ (70 ) 

(G%)m 4a 3 (/? 3 ) m * 3 2 sin(a3>0cos[^ 



2 

3 where 



6 

7 and 
8 



9 CA)-=V(A)i+*2 2 -^ 2 . < 72 > 

10 

11 and the subscript 3 denotes the third experiment. The smallest 

12 residual value corresponds to the correct value of index m and 

13 the correct values of Lame constants. Poisson' s ratio is then 

14 calculated using 
15 

16 o = - 1 . (73) 

2(// + >l) 

17 

18 Young' s modulus can be calculated with 

19 

I.W 4 "' (74) 

21 
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I and the shear modulus can be determined using 

2 

3 G = M . (75) 

4 

5 The above measurement method can be simulated by means of a 

6 numerical example. Soft rubber-like material properties are 

7 used in this simulation. The material has a Young' s modulus E of 

8 [(.Ie8-i2e7) + (5e3/-i3e2/)] N/m 2 where f is frequency in Hz, Poisson's 

9 ratio o is equal to 0.40 (dimensionless) , densitypis equal to 

10 1200 kg/m 3 , and a thickness h of 0.01 m. A compressional 

II (acoustic) wave velocity of Cf of 343 m/s for air is used. All 

12 other parameters can be calculated from these values. The 

13 insonif ication angles of zero, twenty, and forty degrees are 

14 chosen to illustrate this method. 

15 FIGS. 3 and 4 are plots of the transfer functions of 

16 equation (57) at zero (x symbol), twenty (o symbol), and forty 

17 degree (+ symbol) insonif ication angles versus the frequency. 

18 FIG. 3 represents the magnitude of the transfer function versus 

19 the frequency and FIG. 4 represents the phase angle versus the 

20 frequency. 

21 Once the transfer functions are known (typically by 

22 measurement but here by numerical simulation) , the dilatational 

23 wavespeed can be estimated using equations (59) - (62) . FIG. 5 
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1 ' is a plot of the function s versus the frequency. FIGS. 6 and 7 

2 are plots of the actual dilatational wavespeed (solid line) and 

3 the estimated dilatational wavespeed (o symbol) versus the 

4 frequency. FIG. 6 depicts the real component and FIG. 7 depicts 

5 the imaginary component. 

6 FIGS. 8 and 9 are plots of the surface defined in equation 

7 (64) versus real and imaginary components of J3 2 at 1800 Hz. 

8 FIG. 8 depicts a gray scale image of the magnitude versus the 

9 real part of p 2 and FIG - 9 depicts a contour plot of the surface 

10 versus both the real and imaginary parts of P 2 . For both 

11 figures, there are six distinct local minima that are labeled in 

12 bold numbers. The seventh local minima corresponds to /? 2 = 0 

13 which implies there is no shear wave propagation; a physically 

14 unrealizable condition at nonzero wavenumber. These six local 

15 minima are processed at a third measurement location according 

16 to equation (70). The results are listed in Table 1. Local 

17 minimum number 3 has the smallest residual value and corresponds 

18 to the shear wave propagation. The value for ($2)3 is equal to 

19 61.3 + 5. 9i compared to the actual value of fi 2 which is 61.0 + 

20 5.9i. The small difference between the two values can be 

21 attributed to discritization of the surface shown in FIG. 9. 
22 

23 
24 
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Table 1. Values of (fii) m and (e 3 ) m at the Local Minima 



Local Minima 


vaxue oi 




Number m 


\Pl)m 


( hquation 
(70) ) 


1 


22.6 + 5.2i 


0.257 


2 


38.8 + 2.5i 


2.064 


3 


61.3 + 5.9i 


0.013 


4 


94.5 + l.li 


0.426 


5 


125.1 + l.Oi 


0.326 


6 


157.5 + l.li 


0.349 



FIGS. 10 and 11 are plots of the actual shear wavespeed 
(solid line) and the estimated shear wavespeed (o symbol) versus 
the frequency. FIG. 10 depicts the real component and FIG. 11 
depicts the imaginary component. As in FIGS. 8 and 9, the small 
difference between the two values can be attributed to 
discritization of the surface. 

FIGS. 12 and 13 are plots of the actual shear modulus 
(solid line) and the estimated shear modulus (o symbol) versus 
the frequency. FIG. 12 depicts the real component and FIG. 13 
depicts the imaginary component. 
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1 FIGS. 14 and 15 are plots of the actual Young's modulus 

2 (solid line) and the estimated Young's modulus (o symbol) versus 

3 the frequency. FIG. 14 depicts the real component and FIG. 15 

4 depicts the imaginary component. Finally, FIG. 16 is a plot of 

5 the actual Poisson's ratio (solid line) and the estimated 

6 Poisson's ratio (o symbol) versus frequency. Because the 

7 numerical example is formulated using a Poisson's ratio that is 

8 strictly real, no imaginary component is shown in this plot. 

9 Imaginary values of Poisson's ratio are possible and have been 

10 shown to theoretically exist (See T. Pritz, "Frequency 

11 Dependencies of Complex Moduli and Complex Poisson's Ratio or 

12 Real Solid Materials," Journal of Sound and Vibration, Volume 

13 214, Number 1, 1998, pp. 83-104). 

14 The major advantages of this new method is the ability to 

15 estimate complex dilatational and shear wavespeeds of a material 

16 that is slab-shaped and subjected to insonif ication; the ability 

17 to estimate complex Lame constants of the material; the ability 

18 to estimate complex Young's and shear moduli of the material and 

19 the ability to estimate complex Poisson's ratio of the material. 

20 Thus by the present invention its objects and advantages 

21 are realized and although preferred embodiments have been 

22 disclosed and described in detail herein, its scope should be 

23 determined by that of the appended claims. 
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